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^ . Abstract 

Wightman function, the vacuum expectation values of the field squared and the energy- 
momentum tensor are investigated for a scalar field with general curvature coupling param- 
eter in the geometry of a plate in the de Sitter spacetime. Robin boundary condition for 
the field operator is assumed on the plate. The vacuum expectation values are presented as 
'"pi ■ the sum of two terms. The first one corresponds to the geometry of de Sitter bulk without 

JL , boundaries and the second one is induced by the presence of the plate. We show that for non- 

• conformal fields the vacuum energy-momentum tensor is non-diagonal with the off-diagonal 

component corresponding to the energy flux along the direction perpendicular to the plate. 
In dependence of the parameters, this flux can be either positive or negative. The asymptotic 
| behavior of the field squared, vacuum energy density and stresses near the plate and at large 

J> ■ distances is investigated. 

on : 

-xf . 

T-H ■ 

t^! '■ 1 Introduction 

o : 

De Sitter (dS) spacetime is one of the simplest and most interesting spacetimes allowed by general 
relativity. Quantum field theory in this background has been extensively studied during the past 
| two decades. Much of early interest was motivated by the questions related to the quantization 

of fields propagating on curved backgrounds. dS spacetime has a high degree of symmetry 
and numerous physical problems are exactly solvable on this background. The importance of 
this theoretical work increased by the appearance of the inflationary cosmology scenario pp. 
In most inflationary models, an approximately dS spacetime is employed to solve a number of 
problems in standard cosmology. During an inflationary epoch, quantum fluctuations in the 
inflaton field introduce inhomogeneities which play a central role in the generation of cosmic 
structures from inflation. More recently, astronomical observations of high redshift supernovae, 
galaxy clusters and cosmic microwave background [2] indicate that at the present epoch the 
Universe is accelerating and can be well approximated by a world with a positive cosmological 
constant. If the Universe would accelerate indefinitely, the standard cosmology would lead to an 
asymptotic dS universe. Hence, the investigation of physical effects in dS spacetime is important 
for understanding both the early Universe and its future. Another motivation for investigations 
of dS based quantum theories is related to the holographic duality between quantum gravity on 
dS spacetime and a quantum field theory living on boundary identified with the timelike infinity 
of dS spacetime [3] . 
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The Casimir effect is among the most striking macroscopic manifestations of non-trivial 
properties of the quantum vacuum. The presence of reflecting boundaries alters the zero-point 
modes of a quantized field and shifts the vacuum expectation values of quantities, such as the 
energy density and stresses. As a result, vacuum forces arise acting on constraining boundaries. 
The particular features of these forces depend on the nature of the quantum field, the type of 
spacetime manifold, the boundary geometry and the specific boundary conditions imposed on 
the field. Since the original work by Casimir [4] many theoretical and experimental works have 
been done on this problem (see, e.g., [5]-[TT] and references therein). 

The Casimir effect can be viewed as a polarization of vacuum by boundary conditions. The 
interaction of fluctuating quantum fields with background gravitational fields gives rise to an- 
other type of vacuum polarization (see, for instance, jl2l HI])- Here we will study an example 
with both types of vacuum polarization. Namely, we evaluate the vacuum expectation values 
for the field squared and the energy-momentum tensor of a scalar field with general curvature 
coupling parameter induced by a plane boundary on background of (D + l)-dimensional dS 
spacetime. Previously the Casimir effect on background of dS spacetime described in planar 
coordinates was investigated in Refs. p3] for a conformally coupled massless scalar field. In this 
case the problem is conformally related to the corresponding problem in Minkowski spacetime 
and the vacuum characteristics are generated from those for the Minkowski counterpart mul- 
tiplying by the conformal factor. In particular, for the geometry of a single plate the vacuum 
expectation value of the energy-momentum tensor vanishes. In [15] the vacuum expectation 
value of the energy-momentum tensor for a conformally coupled scalar field is investigated in 
dS spacetime with static coordinates in presence of curved branes on which the field obeys 
the Robin boundary conditions with coordinate dependent coefficients (for investigations of the 
Casimir energy in braneworld models with dS branes see Refs. [H]). In these papers the confor- 
mal relation between dS and Rindler spacetimes and the results for the Rindler counterpart |17j 
were used. More recently, the Casimir densities in dS spacetime with toroidally compactified 
spatial dimensions were investigated in |18j . 

The present paper is organized as follows. In the next section the geometry of the problem 
is described and the corresponding positive frequency Wightman function is evaluated for the 
general case of the curvature coupling parameter and for Robin boundary condition on the plane 
boundary. In quantum field theory on curved backgrounds among the important quantities 
describing the local properties of a quantum field and quantum back-reaction effects are the 
expectation values of the field squared and the energy-momentum tensor. In sections [3] and HI 
by using the expression for the Wightman function, we investigate the parts in these expectation 
values induced by the plate. Simple asymptotic formulae are obtained for small and large proper 
distances from the plate with respect to the dS curvature scale. The special cases of Dirichlet 
and Neumann boundary conditions are discussed in section The main results are summarized 
in section El 

2 Wightman function 

For a scalar field with curvature coupling parameter £ the field equation has the form 



where V; is the covariant derivative operator and R is the Ricci scalar for the background 
spacetime. The values of the curvature coupling parameter £ = and £ = £d = (D — 1)/4D, 
with D being the number of spatial dimensions, correspond to the most important special cases 
of minimally and conformally coupled fields. In the present paper the background geometry is 




(1) 
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the (D + l)-dimensional de Sitter spacetime. We write the corresponding line element in planar 
coordinates most appropriate for cosmological applications: 

ds 2 = dt 2 -e 2t ' a J2(dz i ) 2 . (2) 
i=i 

The Ricci scalar and the corresponding cosmological constant are related to the parameter a in 
the scale factor by the formulae 

R = D(D + l)/a 2 , A = D(D - l)/2a 2 . (3) 

In addition to the synchronous time coordinate t it is convenient to introduce the conformal time 
in accordance with r = —ae~ t / a , — oo < r < 0. In terms of this coordinate the line element 
takes the conformally flat form: 

ds 2 = a 2 T- 2 [dT 2 -Y^{dz 1 ) 2 }. (4) 
i=l 

We will assume the presence of a plate located at z D = on which the field obeys the Robin 
boundary condition 

(l + pd z n)<p = 0, z D = 0, (5) 

with a constant coefficient j3. Robin type conditions are an extension of Dirichlet and Neumann 
boundary conditions and appear in a variety of situations, including the considerations of vacuum 
effects for a confined charged scalar field in external fields, spinor and gauge field theories, 
quantum gravity and supergravity. Robin boundary conditions with coefficients related to the 
curvature scale of the background spacetime naturally arise for scalar and fermion bulk fields 
in braneworld models. In this paper we are interested in the boundary induced effects on 
the vacuum expectation values (VEVs) of the field squared and the energy-momentum tensor 
assuming that the field is prepared in the Bunch-Davies vacuum state [19]. These VEVs are 
obtained from the corresponding Wightman function in the coincidence limit of the arguments. 
In addition, the Wightman function determines the response of the particle detector of Unruh- 
DeWitt type, moving through the vacuum under consideration (see, for instance, [12]). 
In order to evaluate the Wightman function we employ the mode-sum formula 

W(x,x>) = (0\cp(xMx')\0) = J2 Va {x)<p*(a/) t (6) 

17 

where {(p a (x), (x)} is a complete set of solutions to the classical field equation satisfying 
the boundary condition ([5]). The eigenfunctions realizing the Bunch-Davies vacuum state and 
obeying the boundary condition have the form 

<P*(x) = C^H^XKri) cos[k D z D + a(^)]e ik - z , a = (k, k D ), (7) 
with the notations rj = |r|, K = ' k 2 + kfj, 

k = (kx, . . . , ko-i), z = (z\...,z D - 1 ). (8) 
The order of the Hankel function (x) in ([7]) is given by the expression 

v = [D 2 /4 - D(D + 1)£ - m 2 a 2 } 1/2 , (9) 
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and the function a(k£>) is denned by the relation 

(10) 



a 2ia(k D ) _ - 1 

i/3k D + 1" 



With this definition the function (|7|) satisfies the boundary condition at z D = 0. Note that we 
also have the formulae 

cos[a(k D )] = = , sin[a(fcu)] = (11) 

V 1 + ^/? 2 ^1 + ^/? 2 

The coefficient in the expression for the eigenfunctions is found from the orthonormal- 
ization condition 

/poo 
dz J dz D y/\i\g 00 [<p a (x)d T( p* a ,(x) - <pl,(x)d T <p a (x)] = 6(k - V)5(k D - k' D ). (12) 

By using the Wronskian relation for the Hankel functions, it can be seen that 

H^iK^HiV'iKr,) - HiV(Kr,)Hp(Kr,) = — . (13) 

In addition, by making use of the definition of the function a(kjj), we can derive the integration 
formula 

/ dz D cos[k D z D + a{k D )]cos[k' D z D + a(k' D )] = -5(k D -k' D ). (14) 
Jo ^ 

Combining the last two relations with (|12|) . for the normalization coefficient we find 

e {v—v*)m/2 
2(2na) D - 1 ' 



Cl = L^n-T (15) 



Substituting the eigenfunctions (|7|) with the normalization coefficient (|15p into the mode-sum 
formula ([U]) for the Wightman function, one finds 

W(x, x') = -. 8 ^n'] D/ n^ I dk e ik '( z - z ') r du cos[uz D + a(u)] cos[uz D ' + a(u)] 
(2Tr) 1J+i a 1J J Jo 

xK^Vk 2 + u 2 e- m/2 )K u ( V 'Vk 2 + u 2 e™ /2 ). (16) 

where we have written the Hankel functions in terms of the modified Bessel function K u (z). The 
formula (|16|) can also be presented in an equivalent form 

W(x,x') = W dS (x,x') + {27r) D +la D-i J dke^ z - z ^ J du cos[u(z D + z D ') + 2a(u)\ 

xK^Vk 2 + u 2 e^ i/2 )K u {r,'^k 2 + u 2 e ni/2 ), (17) 

where 

W d s(x,x') = (2 ^?Xi / dke ' k ' (Z ' Z,) J" du co ^< zD ~ zD '^ 

xK„{ V Vk 2 + u 2 e-™ t2 )K v (ri^/k 2 + u 2 e™ /2 ), (18) 
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is the Wightman function for the dS spacetime without boundaries and the second term on the 
right is induced by the plate at z D = 0. Two-point function (fTHj) is well investigated in the 
literature [TUJ [2U1 EH ESS E3] (see also [T2]) and can be presented in the form 

Wis(x , ,<) = a - r( D /2 + ,)r( D /2-,) T 

' 2( D+3 )/ 2 7r( D+1 )/ 2 (n 2 - 1) ( ^ ' 
where Pv(u) is the associated Legendre function of the first kind and 

u = -l + ^=^ i- — - — . 20 

2r/7/ v ; 

An alternative form is obtained by using the relation between the associated Legendre function 
and the hypergeometric function (see, for example, |24j ) . 

In order to transform the boundary induced part in (|17|) into a more convenient form we 
write 

2cosM^ + z D >) + 2a(u)] = e «(^+^')!^zl + e -^+0!^±l (21) 

Now, in the boundary induced part of (|17|) we rotate the integration contour over u by the angles 
7r/2 and —it/ 2 for the integrals with the first and second terms on the right-hand side of (|2ip . 
respectively. Assuming that (3^0, the Wightman function is presented in the form 

x {^(w) [7_„(T/y) + Uv'y)] + lUvy) + I-v(w)] K v {jfy)} y=y/ ^ i? , (22) 

where I v [z) is the modified Bessel function of the first kind. In deriving this formula we have 
used the relations K v (ze l ) = e T um K v (z) =F Tril v (z). In the case (3 > the first and second 
terms on the right of (|2ip have simple poles at the points if (3 and — i//3 respectively. After 
the rotation of the integration contour these points should be avoided by small semicircles in 
the right half-plane. The integrals along these semicircles give additional contributions which 
should be added to the expression ()22[) where now the integration over u is understood in the 
sense of the principal value. In order to not complicate formulae, in the discussion below we will 
consider the case (3 ^ 0. Note that in (|22p the integration over the angular part of the vector k 
can be done explicitly and the expression for the Wightman function takes the form 

( m '\D/2 roc 

w(x, x >) = W **M + ^(dA^d-i ] Q d yy 

x {K u ( V y) [I- V (rfy) + I v (rfy)] + [I u {r,y) + I-viw)] K v {r,'y)} 

" D-2 J (D-3)/2(k\z - Z'\) e -<z D +z D ') p u + ! 



dk At 



| z _z'|)(f-3)/2 u pu-l 



u=^Jy' 2j +k 2 



(23) 



where J u (z) is the Bessel function of the first kind. 

For Dirichlet and Neumann boundary condition one has (3 = and 1/(3 = 0, respectively, 
and the integral over k in (|23p is evaluated by making use of the formula [25] 



, D _ 2 ^-3)/2(fc|z - z'|) e -yW^+^') [2 n^Kn/^iyX) 



(fc|z-z'|)(^-3)/ 2 T^Tfc 2 " "Vvr y \D/2-i ' ^ 
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with the notation A = \J\l — z'\ 2 + (zP + z D ') 2 . For the corresponding Wightman functions we 
find 

(nn'\ D / 2 n> 1 - D f°° 



2(2^)^/2+1 A D/2 

x {K v ( m ) [i-Mv) + Wv)] + Mm) + i-Aw)] kMv)} , (25) 

where J=D and J=N for Dirichlet and Neumann boundary conditions and 5^ = 1, 6^> = — 1. 

3 VEV of the field squared 

The VEV of the field squared is obtained from the Wightman function taking the coincidence 
limit of the arguments. This limit is divergent and some renormalization procedure is necessary. 
The important point here is that for points away from the plate the divergences are contained 
in the part corresponding to the boundary-free dS spacetime and the boundary induced part is 
finite. As we have already extracted the first part, the renormalization procedure is reduced to 
the renormalization of the boundary-free dS part which is already done in the literature (see, 
for instance, [191 [20l [2~T] ) . As a result the renormalized VEV of the field squared is presented in 
the form 

^ 2 > = (c/^ds + (^ 2 )pi, (26) 
where ((/? 2 }ds is the VEV in dS spacetime without boundaries and the part 



4(4vr 



,-(D+l)/2 



T((D-l)/2)a 



(<P )p1 = w/n iwon D-i / dyy[I v (y)+I- u (y)}K v (y) dxG(x,y), (27) 



Jy 



is induced by the plate. Here and in the discussion below we use the notation 

G(x, y) = (x 2 - y 2 )^" 3 )/ 2 e -^/, Mn l (28) 

lj X j 71 1 

Due to the maximal symmetry of the dS spacetime the VEV (tp 2 )ds does not depend on the 
spacetime point. Note that the boundary induced VEV is a function on the combinations z D /rj 
and (5/rj only. The ratio z D /rj is the proper distance from the plate measured in units of the 
curvature scale a. 

For a conformally coupled massless scalar field we have £ = £d, v = 1/2 and the modified 
Bessel functions in (|27p are expressed in terms of the elementary functions with the result 
[Iu(y) + I-v{y)} K u {y) = 1/y. By making use of the formula 



OO /"OO 



o dyj dx (x 2 - ^- 3 >/ 2 / (x) = ^° D/ £ /2) [ dr r D - 2 f(r), (29) 
the boundary induced part in the VEV of the field squared is presented in the form 

= , ( ± } T\ i r dxx^-\-^°^. 

^ /p (47r)^/ 2 r(D/2) Jo f3x-l y ' 

Formula (|29p is obtained changing the integration variable to u = \/ x 2 — y 2 and introducing 
polar coordinates in the plane (y,u). Note that the result (|30p could also be directly obtained 
by using the fact that in the special case under consideration the problem is conformally related 
to the corresponding problem for a Robin plate in the Minkowski spacetime |26t 127] . From this 
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relation it follows that (v? 2 ) p i = [aC??)] 1 £> (¥' 2 )pi j with the scale factor a(rj) = a/r], which leads 
to the result (|30l) . 

Now let us investigate the behavior of the boundary induced VEV {f 2 ) p \ in the asymptotic 
regions of the ratio z D /r]. For small values of this ratio, z D /rj <C 1, which correspond to 
small proper distances from the plate with respect to the curvature scale a, we introduce new 
integration variables z = xz D jr\ and u = yz D jr\. The argument of the modified Bessel functions 



becomes ur\j ' z D and is large. By using the asymptotic relation [I v (y) + I- U {y)\ K v (y) 
large values y and fixed u, we find 



•pi 



{az D lri) x - D 
(4 7r )D/2 r ( jD /2) 



2;r 



.fix/z D + 1 
fix/ z 



D 



/r] <C 1. 



1/y for 



(31) 



Hence, in the limit under consideration to the leading order the VEV coincides with the cor- 
responding result for a conformally coupled massless scalar. For fixed value z D this limit cor- 
responds to early stages of the cosmological expansion. As the boundary-free part (</? 2 )dS is 
a constant, we conclude that at small proper distances from the plate the total VEV of the 
field squared is dominated by the boundary induced part. From (|31h it follows that for Dirich- 
let boundary condition (/? = 0) this part is negative. For non-Dirichlet boundary conditions, 
assuming that z D < |/?|, near the plate the VEV of the field squared is positive. 

In the opposite limit of large values for the ratio z D /rj (large proper distances from the plate 
with respect to the curvature scale a), the main contribution in the integral over x in (|27p comes 
from the region near the lower limit of the integration and for positive values of the parameter 
v to the leading order one has 



_ 5pT((D + l)/2 - 2v) T{v) 
)pl ~ 2 1 -^(2vr)( D + 1 )/ 2 a^- 1 T(l - i 



2z 



D 



D-2u 



, z D h » 1, 



(32) 



where 8p = 1 for = and 6p = — 1 for ^ 0. In the case v = the VEV of the field 
squared behaves as (v? 2 ) P i oc (r]/z D ) D hi{z D /if). For fixed values of z D this limit corresponds to 
late stages of the cosmological expansion. For imaginary values of u, by using the relation 



K v {y) [I- V {y) + I v {y)\ 



7T 



sinh(|i/|7r 



-Im 



(v/2) 



-2i\v\ 



r 2 (i -%\v\) 



(33) 



for small values of y, in the leading order we find 



<^ 2 >pi 



■sin[2|^| ln(2^/?7) + 



a D - 1 (2z D /t]) d sinh(|i/|7r) 
In this formula the constants Bq and 0o are defined by the relation 

2 2i Hr((D + l)/2-2iM) 



(34) 



B e 



-4a 



r 2 (i 



As we see, in this case the behavior of the boundary induced VEV as a function of the proper 
distance from the plate is oscillatory damping. Note that in terms of the synchronous time 
coordinate the expression in the phase of (|34p is written in the form ]n(2z D /i]) = t / a+ln(2z D / a) . 
From the asymptotic expressions given above we conclude that at large proper distances from 
the plate the VEV of the field squared is dominated by the boundary-free part (</? 2 }dS- 
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4 Energy-momentum tensor 

Now we turn to the investigation of the VEV for the energy-momentum tensor. Having the 
Wightman function we can evaluate this VEV by making use of the formula 



(Tik) = lim did' k W(x, x') + (£ - 1/4) g ik ViV l - £V;V fc - £R lk (ip 2 ), 



(35) 



where Ri k = Dgi k /a 2 is the Ricci tensor for the dS spacetime. Similar to the case of the field 
squared, the VEV of the energy-momentum tensor is presented in the decomposed form 



{Tik) — {Tik) as + {Tik) P b 



(36) 



where {Tik) dS is the corresponding VEV in dS spacetime without boundaries and {Tik) p \ is 
induced by the plate. For points away the plate the latter is finite and the renormalization is 
necessary for the boundary-free part only. Due to the dS invariance of the Bunch-Davies vacuum 
this part is proportional to the metric tensor with a constant coefficient and is well investigated 
in the literature |19] [20] 121] . For this reason in the discussion below we will be focused on the 
boundary induced part. 

In order to evaluate the boundary induced part in the VEV of the energy-momentum tensor 
we need the expression for the covariant d'Alembertian acted on the VEV of the field squared: 



ViVV> pl 



4(4vr 



,-{D+l)/2 



T((D - 1)/2)«d+i Jo 



dyy 



3-D 



dx G(x, y) 



1 - D 



-d„ 



4x z 



V 9 

y y J 



where the notations 



K u (z) = z D / 2 K u (z), I v {z) = z°l 2 [I„(z) + I- U (z)] 



(37) 



(38) 



are introduced. Now, by using the formula (|35p . after long calculations, we present the VEV of 
the energy-momentum tensor in the form 



4(4tt) 



T((D - l)/2)a™ Jo 



dyy 



l-D 



dx G(x, y) F l k {y) + x 2 F l k I v {y)K u {y) 



(39) 



where the function G(x,y) is defined by (f28j) . In formula ([39]) we have introduced the notations 
(no summation over i) 



F D (y) 



y 2 i'MK(y) 



jy 2 d 2 + D(£ - i D )yd y + ZD 



FSiy) + 



D — 1 



I u (y)K u (y), i = l,2,...,D-l, 



(40) 



£(2 — D) + 



D-l 



yd y -ZD\l u {y)K u (y), 



and 



2x[(e-l/4) yay + Z]I v (y)K v (y), 

1 



1 - 4£, Ft = 1 - 4£ 



D - 1 



,f d = fE = o. 



(41) 
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As in the case of the field squared, the VEVs of the energy-momentum tensor are functions of 
the ratios z D /rj and (3/rj only. Note that, by using the well-known properties of the modified 
Bessel functions, the function Fy(y) can also be written in the form 

Fo(v) = [{y 2 /^)d 2 y - D(t + tn)yd v -y 2 + D 2 £ + mV] I u (y)K u {y). (42) 

As we see, the vacuum energy-momentum tensor is non-diagonal with the off-diagonal component 
(Tq 3 ) p i which describes the energy flux along the direction normal to the plate. This flux can 
be either positive or negative (see examples plotted in figures below for Dirichlet boundary 
condition). Note that this type of the energy flux appears also in the geometry of a cosmic 
string on background of Friedmann- Robertson- Walker and dS spacetimes [28|. 129] . 

It can be checked that the boundary induced parts in the VEV of the energy-momentum 
tensor satisfy the trace relation 

(2?>p, = [D(t - ^)V,V l + m 2 ] (ip 2 ) pl . (43) 

In particular, the plate induced part is traceless for a conformally coupled massless scalar field. 
Trace anomaly is contained in the boundary-free part only. 

For a conformally coupled massless field we have v = 1/2. In this case I v (y)K v {y) = y ^ 1 , 
and it can be checked that F D (y) = Fg(y) = 0. For the energy density the corresponding 
expression takes the form 

A(ATv\-( D + l )/ 2 r, D + l r°° r°° , 

{+o) P i - DT{ ( D _ 1)/2)aD+ i J o d y I dx ^ x y) 



y 



Xe X J^l {X DV) ' (44) 

and for the stresses along the directions parallel to the plate we have (no summation over i) 
(7?) P 1 = —(D - l)(T °) p i, i = 1,2,... ,D - 1. Introducing in p]) a new integration variable 
z = \J x 1 — y 2 and polar coordinates y = rcos#, z = rsin#, we can see that the integral over 
9 vanishes. Hence, for a conformally coupled field the plate induced part vanishes. This result 
can also be obtained by using the conformal relation between the problem under consideration 
and the corresponding problem in the Minkowski spacetime. In the latter geometry the Casimir 
energy-momentum tensor for a conformally coupled massless scalar field vanishes. 

The general formulae for the VEV of the energy-momentum tensor are simplified in the 
asymptotic regions of the ratio z D jr\. For small values of this parameter, corresponding to 
small proper distances from the plate, by using the asymptotic formulae for the modified Bessel 
functions for large values of the argument, we find the following asymptotic behavior: 

«'>p> ** - (^B r *'°)>M 

< T ^' - I^M/A^.!^^ (45) 
where we have introduced the notation 

_ . . (4tt)- d / 2 f°° . „ „ 2x zx + l 



T(D/2) J zx-1 

Hence, the VEVs diverge on the boundary. The ratio az D /rj in these formulae is the proper 
distance from the plate and the leading divergences in the components (T^) p i, I = 0, 1, . . . , D— 1, 
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are the same as those in the corresponding problem on the Minkowski bulk. For the latter 
problem the components (T®)^ and (T^) p \ vanish. In the case of the dS bulk near the plate 
we have |(T^) p i| <C K^o^pil *C |(Tq) p i|. Note that for a conformally coupled field the leading 
terms vanish and the next terms in the corresponding asymptotic expansions should be taken 
into account. The limit under consideration corresponds to points near the plate or to late 
stages of the cosmological expansion. In this limit the total VEV of the energy-momentum 
tensor is dominated by the boundary induced part. Note that the function 2^(/3/ ' z D ) is negative 
for Dirichlet boundary condition (/? = 0) and is positive for non-Dirichlet boundary conditions 
assuming that z D <C \/3\. 

Now let us consider the behavior of the vacuum energy-momentum tensor in the asymptotic 
region z D /rf 3> 1. In this limit the main contribution in the integral over x in (j39[) comes from 
the region of the integration near the lower limit. As for the field squared, we consider the cases 
of the real and imaginary values v separately. For real v > 0, by using the asymptotic formulae 
for the modified Bessel functions for small values of the argument, in the leading order we find 
(no summation over i) 

. . r(i + f) 2 2 ^r(p + i)/2-2^ ) 

{ ' /pl ~ r(l - v) (2n)( D + 1 )/ 2 a D + l (2z D /r ] ) D -^ ' 1 J 

for the diagonal components and 

D „ 2^T(u) 5pT{(D + ?,)/2-2v)F { u D) 

for the off-diagonal component with 8/3 defined after formula (|32p . In these expressions the 
notations 

F (o) _ DC\ D ~ 2v ^ 
?+ 4 2v ' 

FP = ^D-2u + l)-^-^,l = l,...,D, (48) 

are introduced. As we see, in the limit under consideration the vacuum stresses are isotropic 
and |(T D ) p i| <C |(Tq) p i|. For v = one has the asymptotic expressions 

o ^ 2Dz D D 2 2 2^r((J) + l)/2) 111(2^/77) 

(T ° )pl * 'WTT) (r ° )pl w m a (2n)iD + m a D +H2z D /ri) D^ ( 49 ) 

and for the diagonal stresses we have (no summation over i) |(T?) p i| <C |(^o')pi|, i = 1, ... ,D. 

For imaginary values v we use the asymptotic relation ()33f) . To the leading order this gives 
the results (no summation over I) 

, ,. i„(2ir)-( D+1 >/ 2 2\v\xBV> r , , , , n , . , , 

w%> - - aro sm t 2 '" 1 ln(2 ' / "> + *] ■ 

«>■>' " asra sm ' 2H M2 * ^ + *">] • (50) 

where the constants are defined by the relations 



B 0O e *, = 2 2,H »r(QP + i)/2-2zH) (o 

- ■ r 2 (i-f|^|) w 

^(00)^ _ aM r((i3 + 3)/2-2i|H) (g) 

6 " r2(i-*H) *h- 



(51) 



In this case the damping of the Casimir densities is oscillatory. For a fixed value of z D the 
oscillating period in terms of the synchronous time coordinate is equal to a/(2ir). 
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5 Dirichlet and Neumann boundary conditions 



The formulae given above for the boundary induced parts in the VEVs of the field squared 
and the energy-momentum tensor are further simplified in the special cases of Dirichlet and 
Neumann boundary conditions. By using the formula 



y 



for the VEV of the field squared from (|27p one finds 



with J=D,N for Dirichlet and Neumann boundary conditions and the definition for 5^ is given 
after formula (1250 . 

In the similar way, for the off-diagonal component of the VEV of the energy-momentum 
tensor we have 

n en A&WcT 13 - 1 f 00 K D/2 (2z D y/ri) 



For the diagonal components we find the expressions (no summation over i 
(J ) _ 2*(-')q-»-' r , f/ ^ 

+Fil v {y)k v {y) 



_if D/2+l(^) K D / 2 {z) 



z D/2-l z D/2 



(55) 

2=22:-° y/r? 



with the notations (|4ip . The corresponding asymptotic expressions are directly obtained from 
those given above for the general Robin case. 

In figures Q] and [2] for conformally and minimally coupled D = 3 Dirichlet scalar fields we 
have plotted the plate induced parts in the VEVs of the energy density and the energy flux along 
the direction normal to the plate as functions of the ratio z D /rj. Recall that the latter is the 
proper distance from the plate measured in units of the dS curvature scale a. The numbers near 
the curves correspond to the values of the parameter ma. We have taken these values in the way 
to have both real and imaginary values for the parameter v. In the first case (graphs for £ = 1/6, 
ma = 1/4 and £ = 0, ma = 3/2) the VEVs decay at large distances monotonically, whereas in 
the second case (graphs for £ = 1/6, ma = 1 and £ = 0, ma = 2) the corresponding behavior 
is damping oscillatory (see formulae (JM|), (|4T|) . ([50]) for the corresponding asymptotics). In the 
case of a conformally coupled field with ma = 1 the first zero is at z D /i] pa 7.19 for the energy 
density and at z D /r) ps 9.64 for the energy flux. For a minimally coupled field with ma = 2 the 
first two zeroes are located at z D /r) ps 2.82, 8.67 for the energy density and at z D /r] ps 3.5, 11.3 
for the energy flux. The energy flux near the plate is positive for both conformally and minimally 
coupled scalars. For Neumann boundary conditions the VEVs have opposite signs. 

It is also of interest to see the dependence of the boundary induced parts in the VEVs on the 
mass of the field. In figure [3] we have presented this dependence for the energy density and flux 
in the case of a conformally coupled field for the fixed distance from the plate corresponding to 
z D jr\ = 3. For large values of the mass the corresponding asymptotic behavior is described by 
formulae (1501). 
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2 4 6 8 2 4 6 8 



Figure 1: The boundary induced part in the VEV of the energy density for D = 3 scalar 
field with Dirichlet boundary condition as a function of the ratio z D /rj. The left/right panel 
corresponds to fields with conformal/minimal coupling. The numbers near the curves are the 
values of the parameter am. 




1 2 3 4 5 2 4 6 8 



Figure 2: The same as in figured] for the energy flux. 
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Figure 3: The boundary induced VEVs of the energy density and flux versus am in the case of 
D = 3 conformally coupled scalar field with Dirichlet boundary condition for the fixed distance 
from the plate z D /rj = 3. 

6 Conclusion 

In the present paper we have investigated the VEVs of the field squared and the energy- 
momentum tensor for a massive scalar field with general curvature coupling parameter induced 
by a plane boundary in dS spacetime. The Robin boundary condition was assumed on the 
boundary. As the first step we have constructed the corresponding positive frequency Wight- 
man function. This function is presented in the form of the sum of boundary-free and boundary 
induced parts. For points away from the boundary the latter is finite in the coincidence limit of 
the arguments and can be directly used for the evaluation of the VEVs of the field squared and 
the energy-momentum tensor. These VEVs are decomposed into boundary-free dS and plate 
induced parts. Due to the maximal symmetry of dS spacetime and the dS invariance of the 
Bunch-Davies vacuum the boundary-free parts do not depend on the spacetime point. These 
parts are well investigated in the literature and we were focused on the plate induced parts. 
The latter are given by formula (|27p for the field squared and by formula (|39|) for the energy- 
momentum tensor and they are functions on the combinations z D ' jr\ and P/rj only. The first 
of these ratios is the proper distance from the plate measured in units of the curvature scale 
a. An interesting feature is that the vacuum energy-momentum tensor is non-diagonal with the 
off-diagonal component corresponding to the energy flux along the direction normal to the plate. 
In dependence of the parameters, this flux can be either positive or negative. For special case of 
a conformally coupled massless scalar field the general formula for the VEV of the field squared 
is simplified to (|30p and the plate induced energy-momentum tensor vanishes. These results can 
also be obtained from the corresponding flat spacetime results by using the conformal relation 
between the geometries. 

For general values of the curvature coupling parameter the formulae for the plate induced 
VEVs are further simplified in the asymptotic regions of small and large proper distances from 
the plate. In the first case the leading term in the asymptotic expansions have the form (I3ip for 
the field squared and the form (|45p for the energy- momentum tensor. Near the plate we have 

\(Td)pi\ < \( t d)p\\ < l( r o°)pil and the total VEV of the 

energy-momentum tensor is dominated 
by the boundary induced part. At large proper distances from the plate the behavior of the 
plate induced parts is qualitatively different for real and pure imaginary values of the parameter 
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v. In the first case these parts decay monotonically as (z D /rj) 2u ~ D for the field squared and 
for the diagonal components of the energy-momentum tensor and like (z D /r]) 2u ~' D ^ 1 for the 
off-diagonal component (T c P) p i. In this limit the diagonal vacuum stresses are isotropic. For 
imaginary values of the parameter v the asymptotic behavior of the plate induced VEVs at 
large proper distances from the plate is described by formulae (|34p and (|50p and is damping 
oscillatory. For a fixed value of z D the oscillating period in terms of the synchronous time 
coordinate is equal to a/(2ir). 

The special cases of Dirichlet and Neumann boundary conditions are considered in section 
[H For these cases the formulae for the expectation values are further simplified to (|53h - (l55p and 
they have opposite signs for Dirichlet and Neumann scalars. 
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